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We study the spectral and transport properties of magnons in a model of a disordered magnet 
called Mattis glass, at vanishing average magnetization. We find that in two dimensional space, 
the magnons are localized with the localization length which diverges as a power of frequency at 
small frequencies. In three dimensional space, the long wavelength magnons are delocalized. In the 
delocalized regime in 3d (and also in 2d in a box whose size is smaller than the relevant localization 
length scale) the magnons move diffusively. The diffusion constant diverges at small frequencies. 
However, the divergence is slow enough so that the thermal conductivity of a Mattis glass is finite, 
and we evaluate it in this paper. This situation can be contrasted with that of phonons in structural 
glasses whose contribution to thermal conductivity is known to diverge (when inelastic scattering is 
neglected) . 

PACS numbers: 



I. INTRODUCTION 

For two reasons, the study of low energy excitations in 
spin glasses - spin glass magnons - is met with notorious 
difficulties: first, the energetic frustration characteris- 
tic for glasses implies the existence of many nearly de- 
generate minima, i.e. the conceptual status of "small" 
fluctuations forming on top of any one of those config- 
urations remains somewhat dubious. Second, even if a 
sufficiently inert extremal configuration was known, the 
solution of the appropriately linearized problem would 
still pose a highly nontrivial problem. 

To get the problem at least partially under control, 
magnons in spin glasses are commonly described in the 
language of mean field theory [113. Within this ap- 
proach one finds that the dispersion of the excitations 
forming on top of a background of vanishing average 
magnetization 



uj[p) cx IpI, 



(1) 



is linear in analogy to the spinwave dispersion of an- 
tiferromagnets. Here, oj and p denote frequency and 
wave vector of the excitations, respectively. Mean field 
theory further predicts that magnons in a glass with 
nonzero average magnetization m have a ferromagnet- 
like branch of excitations. 



>{p) oc / 



(2) 



However, in view of the difficulties alluded to above, 
there is no reason a priori why mean field theory qual- 
ifies to describe spin glasses at all. For example, it has 
been known for more than a decade that mean field 
theory breaks down below the critical dimension dc — 2 
(see Refs. (Yet above the critical dimension 

arguments can be given which support the mean field 
theory results Eqs. (Q) and ^ for the excitation spec- 
trum.) 



Furthermore, an important question which mean field 
theory completely fails to address is the transport prop- 
erties of magnons: while excitations in conventional 
magnetic materials propagate ballistically, the disor- 
der inherent to spin glasses renders the dynamics of 
magnons diffusive and, eventually, leads to mechanisms 
of localization. Clearly, these phenomena cannot be de- 
scribed in terms of a spatially uniform mean field 0] . 

It is the main objective of the present paper to intro- 
duce an alternative approach to the problem which is 
not burdened by these limitations. Developed in close 
analogy to the field theory approach to electron dy- 
namics in disordered solids, the formalism below can 
be employed to address both spectral and localization 
properties of magnons. On the other hand, it has noth- 
ing to say about the first problem mentioned above, 
identification and analysis of reference ground states. 
For this reason, we chose to introduce the approach on 
a prototypical variant of a spin-glass, the so-called Mat- 
tis glass, for which this problem simply does not exist. 

To prepare the definition of the Mattis glass, let us 
recall that spin glasses i8| are usually described by the 
exchange Hamiltonian 



H 



2 X! '^v^i^j^ 
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where i, j refer to nearest neighbor sites on some lattice, 
Jij are random exchange constants, and S" , a = 1, 2, 3 
is a vector representing the spin. (In this paper we 
assume that |S| = 5 ^ 1 is sufficiently large so that 
the spin system can be treated as classical.) The spin 
equations of motion are then given by. 



(4) 



where e°^^'^ is the usual antisymmetric tensor. 
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As mentioned above, one of the major diffieulties 
hampering analytical progress on the problem posed by 
0) is that the ground state(s) of the Hamiltonian are 
not known. For this reason, we consider here a sim- 
plified variant of (O, the so-called Mattis glass. The 
Mattis glass is defined by, 

jka 

i.e. a spin glass Hamiltonian for which the exchange 
constants Jjk ~ ^jS,k factorize. Here, J is a constant pa- 
rameter setting the energy units, j, k are nearest neigh- 
bor sites of a d-dimensional lattice, and S^j are bimodally 
distributed random variables taking value 1 with prob- 
ability p and —1 with probability 1 — p. This factoriza- 
tion entails that a degenerate family of spin configura- 
tions minimizing the Hamiltonian can readily be written 
down: 

^r^"^ =5sign(C.)n", (6) 

where n is a unit vector with arbitrary direction (due to 
rotational invariance of Eq. ||SJ)). Thus the Mattis glass 
is a model with a disordered ground state but without 
the frustration inherent to usual spin systems. To un- 
derstand the physical features of Mattis glass magnons 
we need to explore the equations of motion lin- 
earized around ©. (Under the presumed condition 
S' ^ 1 anharmonic fluctuations, i.e. magnon interac- 
tions can safely be neglected.) We emphasize that this 
part of the analysis is of similar complexity as the study 
of magnons in "real" spin glasses. Indeed, arguments 
can be given to the effect that magnons in real and Mat- 
tis glasses, resp., behave by and large equivalently (see 
Ref. 0). 

The first to study magnon propagation in the Mat- 
tis glass was David Sherrington. In a series of papers 
he found the spectral density of magnons in the 
Mattis glass and their scattering length. Additionally, 
Ref. is an important work studying magnons in a 
ID Mattis glass. It has to be mentioned that any ran- 
dom one dimensional spin chain with nearest neighbor 
interaction is automatically a Mattis glass, and since in 
ID the scattering length coincides with the localization 
length all the transport properties of Mattis magnons 
in ID are known ^ .IQi] . Because of that, in this pa- 
per we will concentrate on dimensionalities higher than 
one. Before turning to the methodological aspects of 
the analysis, let us summarize our main results and re- 
late them to earlier work on similar problems. 

Summary of Results: In this paper we will character- 
ize the behaviour of magnons in terms of their disper- 
sion relation i-u{p), the frequency dependent localization 
length, l{uj), and the so-called density of frequencies 

1 ^ 

n=l 

where Wn are the frequencies of the magnon modes and 
N is the total number of modes in the glass. Borrow- 
ing terminology from the physics of disordered electron 



systems, we will henceforth refer to p{uj) as the density 
of states (DoS). In a way to be formulated precisely be- 
low, the localization length, l{oj), is a measure for the 
exponential decay of magnon modes at frequency to. As 
measurable observables related to these quantities we 
will consider the thermal conductance of the glass, and 
its specific heat. 

We first note (see the discussion in the end of the next 
section) that if the average magnetization of a Mattis 
glass is nonzero, the behavior of its magnons at very low 
frequencies closely resembles that of phonons in a struc- 
tural glass. Since phonons in structural glasses have al- 
ready been discussed in Ref. 0, in this paper we con- 
centrate mostly on the Mattis glass with zero average 
magnetization, or with p = 1/2. In some instances be- 
low where we do consider p ^ 1/2, this will be explicitly 
stated. 

For the dispersion relation of the un-magnetized Mat- 
tis glass we find 

Re uj (X p, Im to (x p"^ in 3(i, 

Re UJ (X — , Im lu oc — t, in 2d. (o\ 

This is consistent with the discussion in Refs. ^ and 
where it is argued that Im oj cx p'^~^ above the critical 
dimension d > dc — 2. This behavior can be contrasted 
with the behavior of phonons in structural glasses whose 
dispersion relation is 

Re oc p, Im (X p'^^^. (9) 

Eq. @ is often interpreted as a manifestation of Raleigh 
scattering. A similar formula also holds for glasses with 
nonzero average magnetization M (but with the substi- 
tution y/uj instead of u). 

In Ref. |l3 self consistent diagrammatic methods 
were applied to derive the DoS 

p{uj) oc UJ |log(w)| in 2d, 
p{ijj) oc Lu^ in id. 

Our field theoretical analysis below will confirm this 
result. Also notice that Eqs. I|1U|I coincide with the 
DoS deduced from mean field theory Eq. ^ (up to a 
logarithmic prefactor in 2d). 

The central result of this paper regards the localiza- 
tion properties of magnons. We find that in 2d magnons 
are localized on the frequency dependent scale 

l{uj) oc CJ^TS^ . (11) 

(In contrast, phonons in structural glasses or, equiva- 
lently, magnons in Mattis glass with nonzero M, are 
subject to a weaker localization mechanism leading to 
an exponentially diverging localization length 'll"!.) At 
length scales below the localization length the magnons 
move diffusively. In the 3d case we find that the dy- 
namics is diffusive no matter how small the frequency, 
i.e. there is no localization. 

In both 2d and Bd, the diffusion constant diverges as 
the frequency of magnons goes to zero: 

^H^-ir- (12) 
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(Compare this with -D(w) oc ui^"^^^ for structural 
glasses 0). Dispersive observables such as the ther- 
mal conductivity are related to the product of diffusion 
constant and density of states, p{uj)D{uj) (weighted by 
the specific heat, C(w), and integrated over frequency.) 
A glance at Eq. shows that this quantity is an 

integrable function at small lo. As an important con- 
sequence, we find that, independent of dimensionality, 
the thermal conductivity of magnons in a Mattis glass 
is finite. Specifically, for 2d, 



/47rA2j2 
\2h I ak-^T-^ 



(13) 



Here, A is a momentum cutoff, reciprocal to the minimal 
wavelength of magnons, and a measures the correlation 
volume of spins (often a oc A""* for short range corre- 
lated spins). Notice that the conductivity is only weakly 
disorder dependent. This is a consequence of the fact 
that the diffusion constant /DoS scale linearly/inversely 
linear with the disorder strength. (Similar behavior is 
observed, e.g., for quasiparticlcs in high temperature 
superconductors 0.) 

In 3(i, the thermal conductivity is given by the fully 
disorder independent result. 



A fc^T 



(14) 



Owing to the absence of localization, the thermal con- 
ductance is obtained from H14|) by a trivial multiplica- 
tion by the system size: Ohm's law. 

In this paper we also discuss the Mattis glass with 
nonzero magnetization. The magnon's low frequency 
behavior is closely related to that of phonons in a glass. 
Thus the density of states is p(uj) cx in agree- 

ment with Eq. (O . The thermal conductivity is infinite 
because of the weak scattering at small frequencies, as 
can be seen from the frequency dependent diffusion con- 
stant D{uj) oc . Finally, at frequencies higher than 
a certain frequency Wc, the magnons' behavior crosses 
over to that of a glass with vanishing magnetization. 
In 2d, u>c — JM /a, where J is the spin-spin interaction 
constant multiplied by the square of the lattice constant 
and a is the correlation volume of spins, see below. In 
Sd, ujc — JM/{aA), where A is the inverse lattice spac- 
ing. 

This concludes our preliminary summary of results. 
The rest of the paper is organized as follows. In the 
next section we derive the magnon equations of motion 
for Mattis glass. We then formulate the field theory 
approach (section IIII|) and use it discuss magnon local- 
ization and transport properties in 2d and 3d (section 

m) 



II. THE EQUATIONS OF MOTION 

The most straightforward way to introduce magnons 
on a technical level is by Holstein-Primakoff transfor- 
mation For those spins Sj with > 0, we write, 
in the harmonic approximation (site indices suppressed 



for notational transparency), 

S"" = S- flta, S+ w aV2S, w a''V2S, (15) 
while for those whose < we write 

S' = -S + fltfl, S+ « a''V2S, S~ w aV2S, (16) 

where a, are Holstein-Primakoff bosons. Substituting 
in Eq. Q and expanding to second order in a, a^, one 
obtains the quadratic Hamiltonian 



N 



J 



jk=l 



where 



n = 



h T 

r h 



(17) 



(18) 



and h and T are N x N symmetric matrices defined in 
the following way 



''J 3 



2d, hjj^^ 



1 



(l+e,e,+^), (19) 



r 



(1 - e.e.+M) 



Here j + p refer to the nearest neighbor sites of the site 
j. As discussed in Ref. Il3|, the magnon equation of 
motion then takes the form 



N 







N 



(20) 



where l^v is the TV-dimensional identity matrix. To 
proceed, it is convenient to perform a unitary rotation 



N 















where the unitary matrix U is defined as 







2 I S + Ijv 



U, 



(21) 



(22) 



and ^jA: = 5jk£,k- As a result, Eq. assumes the 
block-diagonal form 



J fh-T 






h-T 



tjj ~ LO 



s I ^' 



(23) 



where the disorder-independent combination /i — F is 
but the d-dimensional lattice Laplacian ■ (As a side 
remark we mention that the unitary equivalence 7i 
— A > states the positivity of the operator Ti - for a 
bosonic problem, a necessary stability condition.) Thus 
we find 



J 



(24) 
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Here the two signs refer to the upper/lower half of the 
vector ijj. Physically, they correspond to two magnon 
branches of opposite frequency. In the context of Mattis 
glass Eq. (|24|l first appeared in Ref. Q. As a last step, 
we take the continuum limit of this equation. Trading 
the discrete index j for a continuous coordinate x, we 
arrive at [Tal 



and Ref. can be consulted to find out how that was 
done. (For completeness, the main results had been 
summarized in the previous section.) In what follows 
we concentrate on the complementary case of the Mattis 
glass with vanishing magnetization M — (or p = 1/2) , 
with the exception of section IV, where M 7^ will also 
be considered. 



(25) 



where A is the continuum Laplacian and the discrete 
random variable got replaced by a function ^(x), ran- 
domly taking values +1 or —1 with probabilities p and 
1 — p at different points in space. 

As it is hard to deal with bimodally distributed 
variables, we are going to replace £^{x) according to 
^ ^ M + V{x), where M is the average magnetization 
per spin M = 2p — 1 and V{x) is a random Gaussian 
variable with zero mean and correlation 



{V{x)V{y))^a6{x-y). 



(26) 



Remembering the definition of V as a continuous ver- 
sion of the variables ^fe, we deduce that the parameter 
a has the meaning of the correlation volume of spins in 
the Mattis glass. (For uncorrelated spins, a is of the 
order of the elementary lattice volume.) 

In view of universality of extended random systems 
with respect to changes in the microscopic realization 
of the disorder, we trust that taking a continuum limit 
and modelling the randomness according to H26() are 
permissible simplifications. At any rate, the final form 
of the equations of motion on which our further analysis 
will be based reads as 



+ LoM + wV{x) ) V = 0, (27) 



where we have suppressed the ± sign in front of the 
frequency. 

Structurally, Eq. H27|) bears resemblance to the 
Schrodinger equation of a quantum particle in a random 
potential. The crucial difference is, however, that what 
would have been an energy in that equation is now pro- 
portional to magnetization, and what would have been 
the disorder strength is proportional to the frequency 
of magnons. 

Eg. I|27II is also very similar to the main equation of 
Ref. where phonons in structural glasses were stud- 
ied. However, what is M -I- V{x) in Eq. l|77|l was the 
mass density m{x) in Ref. ^Jli ^is such, it was a 
strictly positive quantity. In the Mattis glass, on the 
other hand, M + V{x) (being a continuum analog of 
the bimodal distribution S,j G {—1,1}) can and must 
be negative for some values of x. Yet the methods em- 
ployed in Ref. ^ required to take M -\-V{x) as a ran- 
dom Gaussian variable centered around M = {m{x)). 
Without discussing whether that would invalidate any 
of the conclusions of Ref. ^Jl with respect to the struc- 
tural glasses, we immediately deduce that the proper- 
ties of magnons in Mattis glasses with nonzero average 
magnetization M have effectively been already derived, 



III. FIELD THEORY 



Green functions 



All relevant information about the eigenvalue prob- 
lem defined by (|20|l and (|27|l is contained in the ad- 
vanced and retarded Green functions 



1 



(28) 



where e > is infinitesimal and we have set the av- 
erage magnetization to zero, M = 0. An equivalent 
(see Appendix), but for our purposes more convenient 
representation of the Green function reads as 



uV - 



J 



A ± ze sign lo 



V. 



(29) 



The two quantities we shall focus on in the following 
are 

Ci(c^) = (G+(c^;x,x)), (30) 
C2(c.,l^;|x-y|) = 

= {G+{u + X, y) {uj - r!/2; y, x)) , 

where the averaging (. . .) is over "disorder" V . As 
usual, the average DoS can be represented in terms of 
the Green functions as 



p{lo) = ImGi([x'). 



(31) 



G2(ci.',p) can be used to calculate diffusion constant of 
magnons. At small <C w and small p, 



G2(c^,17;p) 



47rp(cL') 



D(w)p2 



(32) 



where D{uj) is the diffusion constant of magnons at fre- 
quency UJ. 

The thermal conductivity of the system is then given 
by Ej 



du p{uj)C{lu,T)D{lu), 



(33) 



where G is the specific heat of magnons at frequency to, 

{huof 



C{lu,T) 



kT^ (2sinh^) 



(34) 



To compute the disorder average involved in the def- 
inition of the correlation functions Gi_2, we employ the 
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field theoretical formalism of the nonlinear cr-model, an 
approach that has been met with tremendous success in 
the field of disordered fermion physics. In fact, the ap- 
plication of this field theory to "glassy" problems is by 
no means original. In a pioneering work, John, Som- 
polinsky and Stephen llj attacked the related prob- 
lem of phonon localization in structural glasses by the 
same formalism. However, while at first sight the pres- 
ence (phonons) or absence (Mattis glass magnons) of 
the parameter M in H27|) may appear to be of minor 
significance, the opposite is the case. In fact, essential 
elements of the structure of the theory depend on this 
point, which is why very different physical results are 
obtained (cf. discussion above.) 

To prepare the averaging over the disorder, we con- 
sider the supersymmetric generating functional: 



Z\V] 



i:>*cxp(-5[5',*]). 



where the action S is defined as 

n 



'-A 
2 



AR 



V- 



AR 



(35) 



(36) 



Here, ^ = {5'"}, a = 1, . . . , 8 is an eight-component 
vector field, half of whose components are complex vari- 
ables, while the remaining components are anticommut- 
ing (Grassmann variables). Of the commuting (anti- 
commuting) variables, one half refers to the retarded 
sector of the theory, the other half to the advanced 
sector. The simultaneous presence of both sectors is 
indicated by the presence of the Pauli matrix op- 
erating in advanced-retarded space. While our so far 
counting accounts for four different types of variables 
(commuting-anticommuting/advanced-retarded) , a fur- 
ther doubling of the number of integration variables is 
required |15| | by the time reversal invariance of the prob- 
lem (technically, the fact that we are dealing with a 
symmetric operator.) 

Referring for a comprehensive introduction to the 
apparatus of supersymmetry in statistical physics to 
Ref. [Tsf we here merely mention that the rational be- 
hind introducing a supersymmetric structure is that 
Z[V] = 1 is automatically unit-normalized. (The op- 
erator determinants resulting from the integration over 
commuting/anticommuting variables, resp., cancel each 
other.) Not having to worry about determinantal pref- 
actors, the average over the Gaussian distribution of the 
disorder becomes a straightforward operation. Again 
referring for a detailed discussion to Ref. we here 
just note that a Hubbard-Stratonovich decoupling of 
the four-fermion 'interaction' generated by the disorder 
averaging results in 

Z = {Z[V]) = J VReM-Sm, (37) 

where R = {i?"^} is an eight-dimensional matrix field, 
whose action reads as 

S[R] = / str(i?2) + istrln5[i?]-i. (38) 



Here, str is the so-called supertrace and 



(39) 



Eq. IpTfjl represents an (admittedly complicated) rep- 
resentation of unity, Z = 1. To make use of the formal- 
ism, we need to relate our correlation functions C1.2 to 
the generating functional. Indeed, it is straightforward 
to verify that the one point correlation function Ci is 
obtained by computing the expression 

Ci(c.) = *((5i(x)5i*(x))^F(x))^, 

where (. . .)^ denotes the functional average over the 
Gaussian action Ijlffifl and refers to the first compo- 
nent of the commuting sector of the field ^. After the 
disorder average this is equivalent to computing 

Ci(w) =aw(S'Hx)5i*(x)^(x)*(x)) . (40) 

As for the 'two-point' function C2, certain care has to 
be exercised with the sign of the frequency arguments: 
phenomena like diffusion and localization are observed 
in the limit of small frequency difference \n\ w. As- 
suming that the offset frequency > is positive, the 
analytic structure (cf. Eq. (|29|l 'l of the Green function 
is thus determined as indicated in H3U|) . Accordingly, 
the two-point function can be represented as 

C2(a;, n- |x - y|) = {S\^) S'* (y) (y) S'* {^) x 

n 



x5'(x) ( uj+'-^cr^" ] *(x) X 



(41) 



x^(y)(a; + ^a3^«)*(y) 



Here and 5^ stand for advanced and retarded com- 
ponents of the commuting sectors of the vector 'if. 

To actually evaluate the functional expectation values 
(|40|l and 1411) . we need to reduce the (exact) reformula- 
tion of the problem, Eqs. H37|) and H38|) . down to a better 
manageable effective low energy field theory. This re- 
duction is achieved by a gradient expansion around the 
spatially uniform stationary phase configurations of the 
action H38(l . These reference states are determined by 
solution of the equation 



6S[R] 



0=0 



SR{yi) 



= 0^ 



(42) 



^ R = auj 



2tt 



where in the second equality p is momentum. (Owing to 
^ oj, the dependence of the mean field configurations 
on the frequency mismatch Q is weak, which is why f2 
is neglected in 

The simplest solutions to Eq. H42|l acquire the matrix- 
diagonal form 

R = qi+iq2<J3^, 
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where qi^2 are real parameters. The detailed form of 
qi,2 depends on the dimensionality of the problem. To 
avoid reiterations, we therefore first outline the general 
(independent of dimensionality) architecture of the the- 
ory before discussing the cases d= 1,2,3 separately. 

All we will rely upon in the following is that, at suffi- 
ciently small w, (7i is parametrically larger than g2- The 
significance of this relation becomes clear when we note 
that, physically, R plays the role of a self-consistent 
Born (SCBA) self-energy of the problem (in diagram- 
matic language, the self energy calculated neglecting 
diagrams with crossing impurity scattering lines.) 

Indeed, (B^ is structurally equivalent to the familiar 
form of the SCBA self energy equation. The statement 
Qi 3> 92 thus implies that the energy shift acquired by 
impurity induced virtual transitions from one magnon 
mode into others exceeds the energy broadening^ i.e. the 
instability of magnons in a non-translationally invariant 
environment. Again alluding to the formal analogy to 
the Green functions of disordered electrons, qi Ep 
plays the role of the "Fermi energy" of the problem, and 
1 / q2 ~ T is the analog of the inverse fermionic scattering 
time, r. The fact that 91/(72 ~ EpT 3> 1 implies that 
we are working with the analog of a weakly disordered 
system. At the same time, the parameter qi/q2 S> 1 
stabilizes the construction of the low energy field theory 
to be outlined momentarily. 



B. Construction of a low energy action 

The fact that in the limit ^ the action H38|l is 
isotropic in the internal matrix space of the theory im- 
plies that R is but a single element of an entire manifold 
of solutions of the saddle point equation. Indeed, any 
configuration 



TRT- 



represents another solution. Here T g OSp (4|4) is an 
eight-dimensional supermatrix lying in the supergroup 
OSp (4|4), i.e. the maximal group manifold compatible 
with the internal symmetries of the problem jl 5]. Gen- 
eralizing to slowly fiuctuating rotations, we obtain the 
matrix field 

g(x)EEr(x)4«T-i(x)e 

e OSp (4|4)/OSp (2|2) X OSp (2|2) (43) 

as the central degree of freedom of the theory. 

The final step in the construction is to expand the 
action in slow fluctuations of the field Q(x). The fi- 
nite cost of these fluctuations is due to (a) their spatial 
variation and (b) the presence of the so far neglected 
frequency mismatch, 17. Fortunately, the job of actu- 
ally determining the resulting contributions to the ac- 
tion has been done before [l^. AH we need to do is 
carefully identify the parameters of the present problem 
with those of its electronic analog; the algebraic struc- 
ture of the two problems is almost identical. Indeed, 
substituting the 'soft' configurations R = qi+iq2Q into 



the action we obtain 
S[Q] = strln 



+ (Zi + iq2Q + ^af^ {qi + iq2Q) 



(44) 

It is useful to relate this expression to the action of the 
fermionic problem (see Ref. |0|) 



S^[Q] = strln ( -^A + Ep + + ^a^^ ] . (45) 



27' 



Comparison of the two actions leads to the list of iden- 
tifications summarized in the table I below. 

Outgoing from H45|l . an effective low energy action 



Si«[Q] j d'^r str {D^dQdQ + 2inQ(j^^) 



(46) 



can be derived by leading order expansion in the param- 
eters VtT <^ 1 and 1/ L <^ 1, where L is representative of 
the length scales we wish to probe. Eq. H46|) contains 
the bare values of DoS, v, and diffusion constant, Z?Oj of 
the fermionic problem. These quantities are related to 
the parameters of the prototypical action, Eq. H45II as 
also summarized in the table I. Details of the derivation 
can be found in Ref. . 

Quite analogously, starting from the action Eq. (|44|l 
the effective low energy action of the magnon problem 
can be obtained as 



9f 



{2TTYSq2 \ J 



I -9192 AR 



dQdQ+ 
(47) 



We next proceed to discuss what can be learned from 
this representation of the problem. The first quantity 
we would like to calculate is the magnon density of 
states p{uj). This can be extracted from the correla- 
tion function Ci with the help of Eq. H31|) using the 
relations Eq. gUJl and Eq. g^l- The result is 



4(7192 



(48) 



Using Eq. (|48|l . we rewrite Eq. H47|) in a form similar to 
the effective action of the fermion problem, Eq. (|46(l : 



ScB — 



■np{uj) 



d'^r [D{uj)dQdQ + 2iVtQa.^ 



Am 

3 



(49) 



Here D[ijj) is the diffusion constant of magnons at fre- 
quency UJ, 



D{u;) 



{2ttY \ J J 4(7; 



(50) 



Finally, the product p{uj)D{uj), relevant for the calcu- 
lation of the thermal conductivity (cf. Eq. I33|l reads 
as 



p{lo)D{lo) 



r('') ql I 2 

(2^g7 \J 



(51) 



7 



quantity 


fermions 


magnons 


Fermi energy 








scattering time 


T 




(2^2)-^ 


mass 


m 




T— 1 

J 


diffusion constant, Do 
density of states, v 


n 2Eft 

Do dm 


-2 

'~ dm 


Do = ^ 
_ r<'') /2<7i\— d 



TABLE I: Comparison of basic (upper part) and derived (lower part) scales of the problem. In the fifth row, v denotes the 
DoS per volume and F''*' is the volume of the d-dimensional unit sphere. 



Calculating the correlation function C2 using Eq. H41|l . 
we indeed reproduce the expression Eq. 

Finally, we note that the dispersion relation of 
magnons is given by the poles of the Green's function 
G{uj, p) and can be extracted from 



J 



-p =qi+iq2- 



(52) 



Since the functional form of the mean field parame- 
ters qi^2 depend on the dimensionality of space, we next 
consider the cases d = 1,2,3 separately. 



IV. RESULTS 
A. Id case 

Strictly one-dimensional Mattis glasses are amenable 
to transfer matrix techniques and have been discussed 
before us Nevertheless we would like to try to an- 
alyze them here using the techniques proposed in this 
paper. Our conclusion, however, is going to be that, as 
always for strictly one dimensional problems, the sigma 
model description breaks down. 

In d = 1, the momentum integral appearing in the 
saddle point equation (|42ll can readily be done and we 
obtain 



R 



2 4\ 1/3 



2J 



3(7r/3) (1 + i tan(7r/3)cr^^) . (53) 



Substitution of this result into (|^ leads to an SCBA 
DoS 



(54) 



divergent at small frequencies. Eq. H54|l in fact coin- 
cides with the exact expression for the density of states 
obtained in Ref. by the transfer matrix formalism. 
A direct application of SCBA to the ID Mattis glass 
was attempted even earlier, in Ref. fio|. However, from 
a purist point of view, the result does not seem quite 
trustworthy. A glance at (|53|l shows that real (gi) and 
imaginary (52) part of the SCBA self energy are actu- 
ally of the same order. On the other hand, the SCBA 
as such is stabilized by the smallness of the parameter 
92/?! • Put differently, in a one-dimensional environ- 
ment the SCBA, and for that matter the construction 
of our low energy field theory, are ill-founded. For this 
reason, we are not going to discuss the case d = 1 any 
further and turn to d = 2 instead. 



B. 2d case 

In d — 2, the saddle point equation H42|l takes the 
form 



R 



(2^ 



d^p 



Ip^-R 



(55) 



Cutting the logarithmic divergence of this integral by 
introducing a momentum cutoff A, we obtain 



it — log 

2ttJ ^ 



AV 



To logarithmic accuracy, this equation is solved by R - 
<Zi + j'720'^^, where 



/ttAV^ 
91 = TT-f log — 
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auj 
~2T 



Notice that we are now on safe ground inasmuch as 



qi/q2 ^ log I — I > 1, 



(56) 



or in other words, as long as the frequency is small 
enough, 



AJ 



/a 



(57) 



Armed with qi and 52 we can now calculate the DoS, the 
diffusion constant and the thermal conductivity. With 
the help of Eq. (|^ we find 



log 



^AV2 



(58) 



This indeed coincides with the density of states derived 
in Ref. ^3 ■ Looking at the dispersion relation Eq. (|52|l 
with the help of q^ ^ qi we see that 



Recj = 



Imw = 



47r 



\alog(Aj) 



J 

■2^' 



^a^log(f^) 



2^- 
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Next, the diffusion constant is given by 



4auj 



(59) 



Tfiis allows us to calculate the thermal conductivity by 
using Eqs. (|S5|) . and (|^ . Quite remarkably, the 

integral in Eq. is convergent at small frequency 
and we obtain a closed formula for the magnon thermal 
conductivity in 2d 



I2h °® I aBT-^ 



(60) 



As a consequence of Eq. H56|l , this formula works if the 
logarithm in it is large, or T <C KJ jk^fa. 

Finally, we would like to calculate the localization 
length of phonons at frequency lo. To do that, we eval- 
uate the = effective field theory action Eq. (|47|l to 
find 



'5'cff = log 



(fr&lxdQdQ. (61) 



In two dimensions, the supersymmetric cr-model is 
known [m to flow to a disordered phase. The magnon 
localization length is the length scale at which the corre- 
sponding RG group equations renormalize the coupling 
constant of the effective field theory. 



1 

32^ 



log 



2 721 



(62) 



down to values of the order of 1. That gives (compare 
with Ref. _11]) 



or 



cx exp<j ^log 



7rAV2 



(63) 



(64) 



The localization length is divergent as a power law 
of the frequency. This is in contrast to the behavior 
of phonons in structural glasses and, by exten- 
sion, of magnons in Mattis glass with nonzero overall 
magnetization. There the localization length diverges 
much faster with decreasing frequency, as an exponen- 
tial of the inverse frequency square. In other words, 
in these systems, static disorder is less efficient a scat- 
tering agent than in the Mattis glass with vanishing 
magnetization. 



At sufficiently small frequencies lo 

qi q2 OC -j= > 1, 



(66) 



and the effective field theory approach works in this case 
as well. The DoS follows as 



p{uj) = 8 



(aA)2 



(67) 



in agreement with the general mean field theory result 
above the critical dimension d > dc — 2. The dispersion 
relation of magnons is 



Recj = 



Imcj = 



7r2j2 
2aA^ 



7r2j2 



4A V 2aA 



in agreement with the general result Rew oc p, Imw oc 



P 



d-l 



The diffusion constant is given by 



D{lu) = 



J3 



6^2 \/Aa3w2 



(68) 



The thermal conductivity can be found with the help 
of Eq. ^ to give 



A fc2T 



(69) 



In Sd, the sigma model at weak bare coupling EpT ~ 
qi/Q2 ^ 1 is known to flow to a metallic phase where 
the conductance K ~ kL is Ohmic. Accordingly, the 
heat conductance of a sample of linear size L will be 
given by 



K 



KL k^T 

97^"^' 



(70) 



Due to the condition, Eq. (|66|l . the temperature has to 
satisfy 



(71) 



D. Quasi-1 d case 



C. 3d case 

Conceptually, the analysis of the 3d cause parallels 
the discussion of the previous section. We therefore re- 
strict ourselves to a brief statement of the key formulae. 

In analogy to the 2d case, we solve the mean field 
equations Eq. 142|) to 

q;lj2 
91 = ^TT^, 



7r2j 



92 = 



j6A 



27r4 J4 



(65) 



Quasi Id systems are highly anisotropic such that the 
extension in one direction (the "longitudinal" direction) 
is far in excess of the "transverse" extensions. Here we 
will consider such a three-dimensional "wire" made of 
Mattis glass. 

Let us denote the width of the wire, when measured 
in units of lattice spacings, N. Then the actual width of 
the wire is given by N/A, where A, as before, is the in- 
verse lattice spacing. N = 1 corresponds to the strictly 
one-dimensional problem, while here we will consider 
the ^ 1 case. 

As a result of having finite N, the transverse momen- 
tum is quantized in units of (roughly) A/N. The saddle 
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point equation Eq. (|42|l must now involve both integra- 
tion over longitudinal momentum and the summation 
over the transverse ones. Fortunately it is possible to 
go back to three dimensional integration if 



92 » J ( ^ 



(72) 



Once the integration is three dimensional, we can bor- 
row q2 from Eq. H65() to find 



(73) 



a 



This is the condition which the frequency and the num- 
ber of channels must satisfy for the Sd density of states 
be applicable to the quasi-lc? geometry. 

At the same time, for the formalism to work the fre- 
quency cannot be larger than Jy/ A/a, due to Eq. Ht)t)|l . 
We see that as a result, the frequency must lie in the 
interval 



1 

Ni 



J V A 



« 1, 



(74) 



which is always possible at large enough N. 

Once the magnon frequency lies in the interval given 
by Eq. H74|l , the magnons are described by the 3d sigma 
model derived in the previous subsection. However, for 
sufficiently narrow wires and long times the diffusion of 
magnons becomes purely one-dimensional. For that to 
happen, we need to take the frequency f2 to be much 
less then the inverse Thouless time D{lo) / {N / . In 
other words, using Eq. (|68|) we find 



At J3 



Ar2 



(75) 



Our conclusion is that in the long quasi-ld wires the 
magnons whose frequencies satisfy Eq. H74(l diffuse and 
localize purely one-dimensionally for times longer than 
the inverse J7 given in Eq. H75|) above, even though 
their density of states and their diffusion constant as 
a function of frequency lo are given by the 2>d expres- 
sions Eqs. (EI) and (EHJ. 



V. NONVANISHING MAGNETIZATION 
A. Phonon-Magnon crossover 

It is also instructive to see how the magnons behave 
when the average magnetization is nonvanishing M > Q. 
While the properties of very low frequency magnons co- 
incide with those of phonons in structural glasses dis- 
cussed in ,11.] . at somewhat higher frequencies the be- 
havior crosses over to the one indistinguishable from 
that of the zero magnetization magnons. It is this phe- 
nomenon that we would like to study in this section. 

The starting point of the analysis is Eq. H27|l with 
M > 0. It is fairly straightforward to repeat the anal- 
ysis of section IIIII of the paper to find the following 



effective action, the analog of Eq. (|44|) 
"J 



S[Q] = strln 
+ ^(T^^ {Muj + qi+iq2Q) 



A + Muj + qi + iq2Q+ (76) 



Accordingly, the self-consistent Born mean field equa- 
tion Eq. H42(l changes to 



R 



1 



(77) 



2tt J .1 ±v^-R- Muj 
The solution of this equation has, as before, the form 



R = qi+ q2<J 



AR 



(78) 



Finally, the structural form of Green functions and ob- 
servables themselves changes. Specifically, the Green's 
functions now have the form 



uj{V + M) + -A ± id signw 



{V + M) 



(79) 

Comparison with the fcrmionic effective action Eq. (|45|l 
shows that the role of the Fermi energy is now played 
by Mlij -\- qi . Two distinct regimes can now be identi- 
fied: for M <C qi{Lu)/u!, the magnon properties are in- 
distinguishable from those of Mattis glass magnons in 
the absence of a magnetized background. In contrast, 
for M ^ qi{Lu)/uj, the magnons resemble phonons in 
structural glasses. As we shall see in a moment, at suffi- 
ciently low frequencies, the second scenario is effectively 
realized. 

To demonstrate how this works, let us solve the 
SCBA equation Eq. ((ZEl) at d 2. We find 



au^ , / A2 J 



2Mu) 



92 = 



2J 



If 



(80) 



(81) 



the magnons are in the phonon regime. This will hap- 
pen if 



JM 



(82) 



Conversely, at cj ^ JM/a, the magnons are effectively 
in the M = regime studied above in this paper. How- 
ever, the condition qi ^ q2 must also apply in order for 
the formalism developed in this paper to work. This is 
equivalent to w ^ AJ/y^. The consistency condition 
is thus AJ/y/a :s> JM/a or 



M < vA^a. 



(83) 



The average magnetization M < 1, while A^a > 1. We 
see that Eq. always applies. 
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Analogously, at d = 3, we find that the magnetization 
M can be neglected if 



U} > 



MJ 



(84) 



Eq. 0. Similarly, in 3d g2 oc resulting in Imcj cx p^, 
again in agreement with Eq. I^. The diffusion constant 
goes as -D = Do{u>) oc l/w"^/^, leading to the same di- 
vergent behavior of D{oj)p{uj) cx 1/ui as in 2d. 



Its consistency with qi ^ q2 gives M ^ VA-^a. 



VI. CONCLUSIONS 



B. The Phonon regime 

At small enough frequency, the magnons behave in a 
way reminiscent of phonons in structural glasses. Even 
though most of the results for those were obtained in 
Ref. [Tl| we would like to present some of them here, 
in part for completeness, and in part since there are 
still some distinctions between phonons and magnons. 
This mostly stems from the fact that the basic equation 
for phonons in glasses, while practically identical with 
Eq. l|77|l . had in place of u. 

The main feature simplifying the magnetized theory 
at small frequency is that M always dominates the ef- 
fective "Fermi" energy. As a result, it is possible to 
neglect qi in all the calculations. The Green's functions 
can be approximately calculated as 



M 



Lj{V + M) + ± ie signtj 



(85) 



Because of this structure, the problem enjoys a complete 
analogy to the problem of disordered fermions, and the 
results of the treatment of the latter can be directly 
translated into the language of the former. This fact 
has already been successfully exploited in Ref. [lH . 

In particular, it is clear that the DoS will coincide 
with what we called u in the previous sections of the 
paper, 



1 



j2^i\ri~) 



d'^p 6 



J 



Mlo 



(86) 



This coincides with the spectrum of ferrimagnets 0, Q . 

Concentrating on the case of two dimensions d = 2, 
we find that as before q2 = auj'^/{2J), which leads to 
the diffusion constant 



Diu)=DQiu;) 



2J^M 



(87) 



The integral ^ dbj p{uj)D{uj) is logarithmically diver- 
gent, leading to the infinite thermal conductivity of 
magnons at nonzero magnetization. A similar phe- 
nomenon is also observed with phonons in structural 
glasses (although the degree of divergence is different 
there). The physical consequence of this latter phe- 
nomenon is that in real magnets with nonvanishing 
magnetization the thermal conductivity must be domi- 
nated by inelastic scattering processes (completely ne- 
glected in this paper) rather than by elastic collisions. 

We finally note that the 2d dispersion relation reads 
Imtj cx p"^, in agreement with the general formula 



In this paper we studied the spectral and transport 
properties of magnons in a Mattis glass with vanishing 
average magnetization. We find that in 3c? their motion 
is diffusive and that results in finite thermal conduc- 
tivity given by Eq. In 2d the thermal conduc- 
tance is also finite and given by Eq. (|6U|I . This could 
be contrasted with the behavior of phonons in struc- 
tural glasses (and with magnons in a Mattis glass with 
nonzero magnetization) , whose contribution to thermal 
conductance is infinite (when phonon-phonon scattering 
is neglected). 

Of course, one may object that the Mattis glass as 
such represents a somewhat artificial model system. We 
believe, however, that in spite of its artificial definition, 
this system captures much of the physics of magnon 
propagation in generic disordered magnets. Arguments 
to the effect that the results derived in this paper may 
be of wider significance will be presented in a forthcom- 
ing publication^. 
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VIII. APPENDIX 



We would like to show that 



ljV + -A±ieV 



ijjV + —A ± ie signw 



in the limit when e is small. Once this is shown, the 
expression for the Green's function Eq. (|29|l follows im- 
mediately. 

The proof is based on the fact that — -jA is a positive 
definite operator. For simplicity, let us first consider an 
analog of Eq. (|88|l written for real numbers as opposed 
to operators. Take 



1 



UJV - 



d ± ive 



(89) 



where w is a real number and d is a positive real number. 
Take w > 0. Then if d is negative, tov—d is also negative, 
and e can be completely neglected. If v is positive, the 
product ve can be replaced by e, since we are interested 
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in the limit e — s- anyway. In other words, regardless 
of the sign of v, 



1 



1 



ujv — d± ive ujv — d± ie 



, w > 



(90) 



as e is taken to zero. For negative uj, the logic can be 
repeated to result in 



1 



1 



UJV — d ± ive UJV — ie 
These can be combined to give 

1 1 



, uj<0. 



d± 



UJV ~ d ± ie signw 



(92) 



at arbitrary sign of uj. 

The operator generalization of Eq. can be de- 
rived by similar reasoning: take w > as in the text 
above. To show that Eq. H88|l holds, we calculate ma- 
trix elements of 



G-'^ =ujV+^A±ieV 



(93) 



in the basis where ujV + is diagonal. Denoting the 
eigenvalues oi ujV + as Xa and the eigenvalues of 
as —^a (Mo > 0), we find 



uiV + ^A±ieV 



ah 



(94) 

where the matrix elements of -^A in the reference basis 

are written as —UaclJ'cUlij, U is some unitary matrix 
and summation over the index c is implied. In what 
follows, it is crucial that UacfJ'cU^a > or in other 
words, that the diagonal entries of a positive matrix 
have to be positive. 

The eigenvalues are generally nonzero. As we tune 
the parameter w, the eigenvalues can go through zero 



one at a time. If all of Aq are nonzero, e can be taken 
to zero directly in Eq. H94|) . and we obtain 



Gab 



Sab 
Xa' 



(95) 



If one of Xa is zero, more care is needed. The off- 
diagonal matrix elements of G can be found as the ratio 
of the appropriate minors of G~j^ to the determinant of 
(91) G^^. It is easy to see that the determinant of G~^ is 
of the order of e, while the minors are at least of the 
order of e as well (or they could be of higher order in e). 
As a result, the off-diagonal matrix elements of Gab are 
either constant as e goes to zero, or vanishing. How- 
ever, being a constant for a particular value of ut (for 
which one of the eigenvalues Xa vanished) and zero at 
all other values of uj is tantamount to being zero (as in 
limj^o ie/(A + ie) = as a function of A). 

As far as the diagonal matrix elements of G are con- 
cerned, they are also given by the ratio of the minors 
to the determinant of G~^. The main contribution to 
both in the limit e — > is the product of the appropri- 
ate diagonal elements of G^^. It is easy to see that in 
that limit the matrix elements of Gaa are equal to 1/Aa 
except when Xa — 0. For that one vanishing eigenvalue, 



Ga 



UJ 



{UacHcUl,^ 



le 



Since Uac/J-cU^a > Oj and uj > 0, and since e is taken to 
zero, we can replace this by 

G - 

^aa . 

in that limit. So indeed, the coefhcient in front of e 
can simply be put to 1 from the very beginning, as in 
Eq. IHHl) for w > 0. 

Repeating the same argument for w < 0, we arrive at 
Eq. H88|) for arbitrary sign of uj. 
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